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Abstract These lectures are intended to provide a brief pedagogical review of dark 
matter for the newcomer to the subject. We begin with a discussion of the astro- 
physical evidence for dark matter. The standard weakly-interacting massive particle 
(WIMP) scenario — the motivation, particle models, and detection techniques — is 
then reviewed. We provide a brief sampling of some recent variations to the stan- 
dard WIMP scenario as well as some alternatives (axions and sterile neutrinos). 
Exercises are provided for the readerQ 



1 Introduction 



Dark matter is an essential ingredient in a good fraction of the literature on extra- 
galactic astronomy and cosmology. Since dark matter cannot be made of any of the 
usual standard-model particles (as we will discuss below), dark matter is also a cen- 
tral focus of elementary-particle physics. The purpose of this review is to provide a 
pedagogical introduction to the principle astrophysical evidence for dark matter and 
to some of the particle candidates. 

Rather than present a comprehensive survey of the vast and increasingly precise 
measurements of the amount and distribution of dark matter, we will present very 
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simple ("squiggly-line") arguments for the existence of dark matter in clusters and 
galaxies, as well as the arguments for why it is nonbaryonic. The motivation will be 
to provide insight into the evidence and arguments, rather than to summarize results 
from the latest state-of-the-art applications of the techniques. 

Likewise, construction of particle-physics models for dark matter has become a 
huge industry, accelerated quite recently, in particular, with anomalous cosmic -ray 
and diffuse-background results lUEl. Again, we will not attempt to survey these 
recent developments and focus instead primarily on the basic arguments for parti- 
cle dark matter In particular, there has developed in the theoretical literature over 
the past twenty years a "standard" weakly-interacting massive particle (WIMP) sce- 
nario, in which the dark-matter particle is a particle that arises in extensions (e.g., 
supersymmetry U or universal extra dimensions pj) of the standard model that 
are thought by many particle theorists to provide the best prospects for new-physics 
discoveries at the Large Hadron Collider (LHC). We therefore describe this basic 
scenario. More detailed reviews of weakly-interacting massive particles, the main 
subject of this article, can be found in Refs. |l3]|5]|6l. 

After describing the standard WIMP scenario, we provide a brief sampling of 
some ideas for "non-minimal" WIMPs, scenarios in which the WIMP is imbued 
with some additional properties, beyond simply those required to account for dark 
matter. We also briefly discuss some other attractive ideas (axions and sterile neu- 
trinos) for WIMPs. Exercises are provided throughout. 



2 Astrophysical evidence 

It has been well established since the 1930s that there is much matter in the Uni- 
verse that is not seen. It has also been long realized, and particularly since the early 
1970s, that much of this matter must be nonbaryonic. The evidence for a significant 
quantity of dark matter accrued from galactic dynamics, the dynamics of galaxy 
clusters, and applications of the cosmic virial theorem. The evidence that much of 
this matter is nonbaryonic came from the discrepancy between the total matter den- 
sity Q.,n ~ 0.2 — 0.3 (in units of the critical density = 2>Hq/%kG, where Hq is 
the Hubble parameter), obtained from such measurements, and the baryon density 
Q,y ~ 0.05 required for the concordance between the observed light-element (H, D, 
^He, ^He, ^Li) abundances with those predicted by big-bang nucleosynthesis IQ, 
the theory for the assembly of light elements in the first minutes after the big bang. 

Rather than review the historical record, we discuss the most compelling argu- 
ments for nonbaryonic dark matter today as well as some observations most relevant 
to astrophysical phenomenology of dark matter today. 
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2.1 Galactic rotation curves 

The flatness of galactic rotation curves has provided evidence for dark matter since 
the 1970's. These measurements are particularly important now not only for estab- 
lishing the existence of dark matter, but particularly for fixing the local dark-matter 
density, relevant for direct detection of dark matter. We live in a typical spiral galaxy, 
the Milky Way, at a distance 8.5 kpc from its center. The visible stars and gas in 
the Milky Way extend out to a distance of about 10 kpc. From the rotation curve, 
the rotational velocity Vc{r) of stars and gas as a function of Galactocentric radius r, 
we can infer the mass M< (r) of the Galaxy enclosed within a radius r. If the visible 
stars and gas provided all the mass in the Galaxy, one would expect that the rotation 
curve should decline at radii larger than the 10 kpc extent of the stellar disk accord- 
ing to the Keplerian relation ^ GM^ts/''- Instead, one observes that \'c{r) remains 
constant (a flat rotation curve) out to much larger radii, indicating that M< (r) °^ r for 
r S> 10 kpc and thus that the Galaxy must contain far more matter than contributed 
by the stars and the gas. 

Assuming a spherically symmetric distribution of matter, the mass inside a radius 
r is given by 
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Fig. 1: Measured rotation curve of NGC6503 with best fit and contributions from halo, disk and 
gas. From Ref (8) 
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M<(r) =4;r / p{r'y^dr' . (1) 





An estimate for the distribution of dark matter in the Galaxy can be obtained from 
the behavior of the rotation curve in the inner and outer galaxy. For example, the 
density distribution for the cored isothermal sphere, given by, 

p{r)=pQ 2.2 ' (2) 
n + 

where /? ~ 8.5 kpc is our distance from the Galactic center and po is the local dark- 
matter density, provides a qualitatively consistent description of the data. For large 
r, p ~ r^^ =^ M{r) oc r =J> v ^ const, while for small r, p ^ const => M{r) ^ 
V cx: r. Eq. ^ describes a 2-parameter family of density profiles and by 
fitting the observed data one finds a scale radius a ^ ?> ~ 5 kpc and local mat- 
ter density po ~ 0.4GeVcm^'^; the uncertainties arise from standard error in the 
rotation-curve measurements and from uncertainties in the contribution of the stel- 
lar disk to the local rotation curve. Because the dark matter is moving in the same 
potential well, the velocity dispersion of the dark matter can be estimated to be 

1 /2 

(^'dm) ^ km/sec. The simplest assumption is that the dark matter has a 
Maxwell-Boltzmann distribution with f{v) ~ e^'"/^*'", where v ~ 220km/sec. 
Exercise 1. Explain/estimate how po would be affected if 

• (a) the halo were fattened, keeping the rotation curve unaltered; 

• (b) the profile were of the Navarro-Frenk-White (NEW) type: p (r) oc p^/ [r{r + 
rc)^], keeping the local rotation speed the same; 

• (c) the stellar contribution to the rotation curve was either increased or decreased. 



2.2 Galaxy Clusters 

Galaxy clusters are the largest gravitationally bound objects in the Universe. They 
were first observed as concentrations of thousands of individual galaxies, and early 
application of the virial theorem ~ GM/R (relating the observed velocity disper- 
sion to the observed radius R of the cluster) suggested that there is more matter in 
clusters than the stellar component can provide i9|. It was later observed that these 
galaxies are embedded in hot x-ray-emitting gas, and we now know that clusters are 
the brightest objects in the x-ray sky. The x rays are produced by hot gas excited to 
virial temperatures T ^ keV of the gravitational potential well of the dark matter, 
galaxies, and gas. A virial temperature T ^ keV corresponds to a typical velocity 
for the galaxies of v ~ 10^ km/s. 

Observations of clusters come from optical and x-ray telescopes and more re- 
cently via the Sunyaev-Zeldovich effect [lOJ. Several independent lines of evidence 
from clusters indicate that the total mass required to explain observations is much 
larger than can be inferred by the observed baryonic content of galaxies and gas. 
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2.2.1 Lensing 

Galaxy clusters exhibit the phenomenon of gravitational lensing lfm[T2l . Because 
the gravitational field of the cluster curves the space around it, light rays emitted 
from objects behind the cluster travel along curved rather than straight paths on their 
way to our telescopes lfT3l . If the lensing is strong enough, there are multiple paths 
from the same object, past the cluster, that arrive at our location in the Universe; this 
results in multiple images of the same object (e.g., a background galaxy or active 
galactic nucleus). Furthermore, because the light from different sides of the same 
galaxy travels along slightly different paths, the images of strongly lensed sources 
are distorted into arcs. For instance, HST observations of Abell 2218 show arcs and 
multiple images as shown in Fig.|2] If the lensing is weak, the images may become 
slightly elongated, even if they are not multiply imaged. 

For a lensing cluster with total mass M and impact parameter d the deflection 
angle is of order 

1/2 



fGM\ 



(3) 



Thus, from measurements of the deflection angle and impact parameter (which can 
be inferred by knowing the redshift to the lensing cluster and source), one can infer 
that the total mass M of a cluster is much larger than the observed baryonic mass 
Mb. 

Exercise 2. Suppose a massive particle with velocity v is incident, with impact pa- 
rameter b, on a fixed deflector of mass M. Calculate the deflection angle (using clas- 
sical physics) due to scattering of this particle via gravitational interaction with the 
deflector. Show that you recover a = (GM/dc^^ in the limit v ^ c, the velocity 
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at which light rays propagate. Actually, the correct general-relativistic calculation 
recovers this expression, but with an extra factor of 2. 

Exercise 3. Estimate the deflection angle a forlensing by a cluster ofM ^ lO^^M© 
and for an impact parameter of 1 Mpc. 



2.2.2 Hydrostatic equilibrium 

In a relaxed cluster, the temperature profile T{r) of gas, as a function of radius r, 
can be inferred using the strength of the emission lines, and the electron number 
density ne(r) can be inferred using the the x-ray luminosity L{r). Combined, these 
observations give an estimate of the radial pressure profile p{r) oc ne{r)kBT{r). In 
steady state, a gravitating gas will satisfy the equation of hydrostatic equilibrium, 

dp ^ M<(r)pgas(r) 

-1- - -<J 2 ■ 

ar n 

Here, M< (r) is the total (dark matter and baryonic gas) mass enclosed by a radius r 
and Pgas('") is the density at radius r. Eq. (|4]i can be used to determine the total mass 
M of the cluster Comparison with the observed baryonic mass Mt, again shows 
that M ^ Mj,. In particular, observations using the x-ray satellites XMM-Newton 
and Chandra indicate that the ratio of baryonic matter to dark matter in clusters 
is Qt>/£2i„ ^1/6. Additional constraints to the cluster-gas distribution can be ob- 
tained from the Sunyaev-Zeldovich (SZ) effect. This is the upscattering of cosmic 
microwave background (CMB) photons by hot electron gas in the cluster; the mag- 
nitude of the observed CMB -temperature change is then proportional to the integral 
of the electron pressure through the cluster (see, e.g., ifTOl ). 

Exercise 4. Estimate, in order of magnitude, the x-ray luminosity Lx for a cluster 
with total mass M ^ lO'^Mj, and a baryon fraction 1/6 in hydrostatic equilibrium 
with maximum radius R ~ Mpc. 

Exercise 5. Assume the cluster in Exercise 4. is isothermal (T{r) = T = const.) 
with a dark-matter distribution consistent with an NEW profile with ~ R/IQ. 
Neglecting the self-gravity of the gas: 

• (a) Show the properly normalized dark-matter density profile is approximately 
p{r) ~ (233/45)M,./[r(r + r,.)2], where Mc=M<{rc) is the mass enclosed within 
the scale radius r^. Determine M< (r) and and in terms ofM for this cluster. 

• (b) Using your results from (a) solve Eq.^and show that the gas density profile 
in such an NEW cluster takes the form Pgas{r) °= (1 + r /rcY'''^l^ , where F <^ 
{GMcHmp/r,)/{kBT). 
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2.2.3 Dynamics 

According to the virial theorem, the velocity dispersion of galaxies is approximately 
(r) ~ GM< {r)/r, where M< (r ) is the mass enclosed within a radius r. Therefore, 
from measurements of the velocity dispersion and size of a cluster (which can be 
determined if the redshift and angular size of the cluster are known), one can infer 
the total mass M. Once again, the total mass is much larger than the baryonic mass 

Cluster measurements are by now well established, with many well-studied and 
very well-modeled clusters, and there is a good agreement of estimates of M from 
dynamics, lensing. X-ray measurements, and the SZ effect. The current state of the 
art actually goes much further: one can now not only establish the existence of dark 
matter, but also map its detailed distribution within the cluster 

Exercise 6. Following Zwicky 03t, use the virial theorem to find an approximate 
formula relating the average mass of a galaxy in a galaxy cluster to the observed size 
and velocity dispersion of the cluster assuming that the system is self- gravitating 
(and assuming only that the observed galaxies contribute to the mass of the system). 
What answer would Zwicky have found for the Coma cluster with modern data? 



2.3 Cosmic Microwave Background and Large-Scale Structure 

Measurements of the cosmic microwave background (CMB) radiation and large- 
scale structure (LSS) of the Universe provide perhaps the most compelling evidence 
that the dark matter is non-baryonic and the most precise measurements of its abun- 
dance. 

One obtains from CMB maps the angular power spectrum Cf of CMB tempera- 
ture anisotropics as a function of multipole i. If the temperature T{h) is measured 
as a function of position n on the sky, then one can obtain the spherical-harmonic 
coefficients a(,„ = f dhT{h)Y^^^^{n). The Q's are then simply the variance of the 
spherical-harmonic coefficients: Q = (\a(i„\^). Theoretical predictions for the power 
spectrum depend on the values of cosmological parameters like the matter density 
Qmh^, the baryon density £2t,h^, the cosmological constant A, the scalar spectral 
index n,, the optical depth T due to reionization, and the Hubble parameter Hq. One 
can thus determine these cosmological parameters by fitting precise measurements 
of the CfS to the theoretical predictions |fT6| . Current measurements provide de- 
tailed information on Cf over the range 2 < I < ^^(lOOO), thus providing precise 
constraints to the cosmological parameters. 

In the year 2000, data from the Boomerang and MAXIMA experiments (with 
supernova measurements) gave i2,„/r =0.13 ±0.05 with error bars that shrink to 
±0.01 taking into account other measurements or assumptions (e.g., LSS, Hubble- 
constant, and supernova measurements, and/or the assumption of a flat Universe) 
IflTll . Now, with WMAP, Q,„h^ = 0.133 ±0.006 and i2fo/r = 0.0227 ±0.0006 [[ISj. 
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Fig. 3: Dependence of the CMB power spectrum on the cosmological parameters. From Ref. 1151 . 



Exercise 7. Suppose that the temperature is measured with a Gaussian noise Ot — 

1 /2 

25 /iK in Npix ^ 10^ pixels on the sky. Estimate the rms temperature (^(5T /T)^'^ 
that results. 



3 Basic properties of dark matter 

Having established the existence of dark matter and presented the case that it is non- 
baryonic, we now consider the requirements for a dark-matter candidate and discuss 
some possibiHties. Every dark-matter candidate should satisfy several requirements: 
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• Dark matter must be dark, in the sense that it must generically have no (or ex- 
tremely weak) interactions with photons; otherwise it might contribute to the 
dimming of quasars, create absorption Unes in the spectra of distant quasars |fT9ll . 
or emit photons. One way to quantify this is by assuming that dark-matter parti- 
cles have a tiny charge fe (where e is the electron charge and / <C 1), which can 
be quantitatively constrained 1201 . 

• Self-interactions of the dark matter should be small. We can estimate the cross 
section for DM-DM scattering in the following way: if DM particles scatter less 
than once in the history of the Universe, then the mean free path is less than 
X=vdmHq^ - (3 X lOW/sec) (lO'^sec) - 3 x lO^^cm. Then, if the galactic- 
halo density is poM ^ lO^^^g/cm^, the opacity for self-scattering in the galactic 
halo is JC = (PdM'^) * = cr/m ^ cm^/g. Thus, if the elastic-scattering cross sec- 
tion is (7 > 10^^^ (OT/GeV)cm^, then K" > 1 and the typical halo-dark-matter 
particle scatters more than once during the history of the Universe. If dark matter 
self-scattered, it would suffer gravothermal catastrophe: that is, in binary inter- 
actions of two dark-matter particles, one particle can get ejected from the halo, 
while the other moves to a lower-energy state at smaller radius. As this occurs 
repeatedly much of the halo evaporates and the remaining halo shrinks. Although 
a variety of arguments can constrain dark-matter self-interactions, stringent and 
very transparent constraints come from observations of the Bullet Cluster, the 
merger of two galaxy clusters, in which it is seen (from gravitational-lensing 
maps of the projected matter density) that the two dark-matter halos have passed 
through each other while the baryonic gas has shocked and is located between 
the two halos 1211 . 

• Interactions with baryons must also be weak. Suppose baryons and dark matter 
interact. As an overdense region collapses to form a galaxy, baryons and dark 
matter would fall together, with photons radiated from this baryon-DM fluid. 
This would result in a baryon-DM disk, in contradiction with the more diffuse 
and extended dark-matter halos that are observed. If DM interacted with baryons 
other than gravitationally in the early Universe, the baryon-photon fluid would 
be effectively heavier (have a higher mass loading relative to radiation pressure) 
even before recombination, so that the baryon acoustic oscillations in the matter 
power spectrum and the CMB angular power spectrum would be modified |22l . 

• Dark matter cannot be made up of Standard Model (SM) particles, since most 
leptons and baryons are charged. The only potentially suitable SM candidate 
is the neutrino, but it cannot be dark matter because of the celebrated Gunn- 
Tremaine bound l23l . which imposes a lower bound on the masses of dark- 
matter particles that decoupled when relativistic. The argument is the following: 
The momentum distribution in the Galactic halo is roughly Maxwell-Boltzmann 
with a momentum uncertainty Ap ^ my{v) ((v) ^ 300 km/sec), while the mean 
spacing between neutrinos is Ax^Uy^^^ ^ {py /niv)^^^^ . The Heisenberg uncer- 
tainty principle gives AxAp> Ti, which translates into a lower bound niy > 50 eV. 
(This Heisenberg bound can actually be improved by a factor of 2 by using ar- 
guments involving conservation of phase space.) Stronger bounds (niy > 300eV) 
can be obtained from dwarf galaxies which have higher phase-space densities. 
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Fig. 4: An example of a Feynman diagram for annihilation of two WIMFs x (neutralinos in this 
case) to fermion-antifermion pairs (where the fermions are either quarks q or leptons /) via ex- 
change of an intermediate-state squark q or slepton I. 



As discussed below, there will be a cosmological density of neutrinos left over 
from the big bang, with a density Qyh^ ~ 0.1 (my/ 10 eV). The neutrinos of mass 
ffiy > 300 eV consistent with the Gunn-Tremaine bound would overdose the Uni- 
verse. Thus, neutrinos are unable to account for the dark matter 



4 Weakly Interacting Massive Particles (WIMPs) 

Perhaps the most attractive dark-matter candidates to have been considered are 
weakly-interacting massive particles (WIMPs). Many theories for new physics at 
the electroweak scale (e.g., supersymmetry, universal extra dimensions) introduce a 
new stable, weakly-interacting particle, with a mass of order M;,, ~ lOOGeV. 
For example, in supersymmetric (SUSY) theories, the WIMP is the neutralino 



a linear combination of the supersymmetric partners of the photon, Z boson, and 
neutral Higgs bosons. Neutralinos are neutral spin-1 /2 Majorana fermions. In the- 
ories with universal extra dimensions there are Kaluza-Klein (KK) states Ykk, 2%k^ 
H^[^, which are neutral KK bosons. The candidates are stable (or quasi-stable; i.e., 
lifetimes greater than the age of the Universe T ^ fj/) and particle-theory models 
suggest masses ~ 10 — lO-'GeV. 

In typical theories two WIMPs can annihilate to SM particles. For example, for 
a neutralino we have the tree-level diagram in Fig.H] where m- j ~ 100 GeV, so that 



hi 



(5) 



CT - a^mZ^Ml - 10-**GeV-2. 
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4.1 WIMP Freezeout in Early Universe 



We now estimate the relic abundance of WIMPs in the standard scenario of thermal 
production (see, e.g., Ref. Il24ll ). In the early Universe, at temperatures T S> Myr, 
WIMPs are in thermal equilibrium and are nearly abundant as lighter particles, like 
photons, quarks, leptons, etc. Their equilibrium abundance is maintained via rapid 
interconversion of %% pairs and particle-antiparticle pairs of Standard Model parti- 
cles. When the temperature falls below the WIMP mass, however, the WIMP abun- 
dances become Boltzmann suppressed, and WIMPs can no longer find each other to 
annihilate. The remaining WIMPs constitute a primordial relic population that still 
exists today. 

We now step through a rough calculation. To do so, we assume that the WIMP is a 
Majorana particle, its own antiparticle (as is the case for the neutralino, for example), 
although the calculation is easily generalized for WIMPs with antiparticles (e.g., KK 
WIMPs). 

The annihilation rate for WIMPs is r{xx ^ ■ ■ ■) ~ where a is 

the cross section for annihilation of two WIMPs to all lighter standard-model parti- 
cles, V is the relative velocity, and the angle brackets denote a thermal average. The 
expansion rate of the Universe is H ^ {SJtGp /3)'''^ ~ T^/Mpi during the radiation 
era, where p °^ T^. In the spirit of "squiggly lines" we have neglected factors like the 
effective number of relativistic degrees of freedom g^, in the expansion rate, which 
the careful reader can restore for a more refined estimate. 

By comparing these two rates, one can identify two different regimes: 

• At early times, when T ^ M^, n-^ °^ and F ^ H: particles scatter and annihi- 
late many times during an Hubble time and this maintains chemical equilibrium. 

• At late times, when T <C M^, <^ j^ftg^'^'X-l^ (note that the chemical potential 
/Zx = in the case of Majorana particles such as the neutralino) and F <^H: there 
can be no annihilations, and the WIMP abundance freezes out (the comoving 
number density becomes constant). 

This sequence of events is illustrated in Fig.|5] which shows the comoving number 
density of WIMPs as a function of the inverse temperature in equilibrium (solid 
curve) and including freezeout (dashed curves). 

Freezeout occurs roughly when F{Tf) ^ H{Tf). For nonrelativistic particles, 

~8x {M^T /iTiy'^^ e^^^l'^ , so the freezeout condition becomes 



(MyTffl^ e-W ~ ^ ^ ^ ~ In 



Mpi 



MpiM]/^ {av) 



(6) 



where the mass parameters are in GeV. Taking (cv) ~ a? /M^, and taking as a first 
guess Tf ^ M^, we finally find 
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where the numerical values are characteristic electroweak-scale parameters (i.e. CT ~ 
10-^GeV"^ Mjf^ - lOOGeV). 
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Fig. 5: Equilibrium (solid curve) and relic abundance (dashed curves) of WIMP particles. From 
Ref. [3. 



At freezeout, the abundance relative to photons is 
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Today we know that 
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(8) 



(9) 



with no explicit dependence on the particle mass. 

We thus obtain the observed abundance ^x^~ ~ 0.1 for a ^ 10"* (0.1 x 10^^ x 
lO-^^GeV^^ which turns out to be nearly exact, even though we 



10"*^)"'GeV"2 _ in-8r.~->'-2 
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have been a bit sloppy. A more precise calculation (including all the factors we have 
dropped) gives 



a remarkable result, as it implies that if there is a new stable particle at the elec- 
troweak scale, it is the dark matter 

As an aside, note that partial-wave unitarity of annihilation cross sections re- 
quires <7 < M^^, which means 

Z (M;^/300TeV) . This thus requires Q^h^ < 
0.1, < lOOTeV, without knowing anything about particle physics f25l. More 
precisely, this bound applies for point particles and does not apply if dark matter 
particles are bound states or solitons. If the interactions are strong, a ~ 1, the bound 
is already saturated. 

Although our arguments have been rough, one finds in SUSY and KK models 
that there are many combinations of reasonable values for the the SUSY or KK 
parameters that provide a WIMP with Q^li^ ^OA for lOGeV <Mx<l TeV. 

Exercise 8. Eq. ITOl ) was derived assuming that the annihilation cross section (av) 
is temperature-independent. Redo the estimate forQ^h^ assuming that (c7v) °^ T" , 
where n~ 1,2,3, ••• . 



4.2 Direct detection 

If WIMPs make up the halo of the Milky Way, then they have a local spatial den- 
sity -0.004 (M;^/100GeV)"'cm-3 (roughly one per liter), and are moving with 
velocities v ~ 200 km sec^'. Moreover, there is a crossing symmetry between 
the annihilation qq and the elastic scattering ^ Xl processes — apart 

from some kinematic factors the diagrams are more or less the same (as shown 
in Fig. |6ll — so one expects roughly that the cross section o{x<l ^ XO) ^ '^{XX ^ 
qq) ~ 10^^^ cm". One can therefore hope to detect a WIMP directly by observing its 
interaction with some target nucleus in a low -background detector composed, e.g., 
of germanium, xenon, silicon, sodium, iodine, or some other element. 




(10) 



X 



q 




X 



q 



Fig. 6: Crossing symmetiy between annihilation and scattering diagrams. 
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At low energies, quarks are bound into nucleons and nucleons in turn are bound 
into nuclei, so the cross section one actually needs is o{xN —> xN) (where here 
stands for a nucleus). The calculation relating the xq interaction to the xN interac- 
tion requires both QCD and nuclear physics. It is complicated but straightforward. 
Here we will simply assume, for illustration, that o{xN X^) ^ '^{XQ ^ X^)- 

The rate at which a nucleus in the detector is hit by halo WIMPs is then 

7?-«;^(7v- (0.004cm"^)(10"^^cm2) (3 x 10^^) - 10"-V"'; (H) 

if there are 6 x lO^^M/{Ag) nuclei in a detector, for an atomic number A ~ 100 we 
expect to see R ^ 10/kg/ yr events. 

Let us estimate the recoil energy of a nucleus struck by a WIMP. If a WIMP 
of ~ lOOGeV runs into a nucleus with A ^ 100, the momentum change is 
4/7 ~ M^v, and the nucleus recoils with an energy of order E ^ {Ap)^/2m ^ 
(100GeV10"^)2(100GeV)-' - lOOkeV. 

To do things more carefully, one has to account for the fact that the cross sec- 
tion one actually needs are the interaction cross sections with nuclei, and via the 
following steps, 

^ix(i)--^(yixn),(y{xp) ,— t . (y{xN), 

QCD nuclear physics 

some theoretical uncertainties are introduced. One also finds that (7{xN) is reduced 
relative to o{xq) by several orders of magnitude. 

Qualitatively, there are two different types of interactions, axial and scalar (or 
spin-dependent and spin-independent). The first is described by the Lagrangian, 

Ui^i^x^^x'fysxqJiiysq. (i2) 

which couples x to the spin of unpaired nucleons; this works only for nuclei with 
spin, and the coupling is different for unpaired protons or neutrons. Through this 
interaction one expects a °^ f', where s is the average spin ^ 1/2 of the unpaired 
proton or neutron in nuclei with odd atomic number 
The second interaction is described by the Lagrangian, 

^scalar °= ZZ??7 (13) 

which couples x to the mass of the nucleus, thus giving a cross section a °^ A^ 
(where M and A are the nuclear mass and atomic number), which implies higher 
cross sections for larger A. However, this scaling is only valid up to a limit. In fact, 
the momentum exchanged is Ap ^ (100GeV)(10^^) ^ 0.1 GeV, and the nuclear 
radius is roue hlyr-A'/^lO-'-^ cm, so from the uncertainty principle one has rAp> 
1 when 

(0.1GeV)(10-'^cm) 
2xlO-i4GeVcm ^ ^ ^ ^ ^^^^ 
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Detailed calculations show that the cross section for WlMP-nucleus elastic scatter- 
ing does not increase much past A > 100. 

In experiments, people usually draw exclusion curves for the WlMP-nucleon 
cross section versus the WIMP mass M^. The exclusion curves are less constraining 
both for low Mj^ because of the low recoil energy, and for large because (for fixed 
local energy density p^j^) the number density oc MZ^ . To date, only the DAMA 
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Fig. 7: Exclusion plot for the spin-independent dark-matter parameter space. The region favored 
by the DAMA annual modulation is inconsistent with the current bound (solid curve) from CDMS. 
The broken curves are forecasts for future experiments. We also show, for illustrative puiposes 
only, predictions for a WIMP model with a lightest-Higgs-boson mass of mi, = 150 GeV. 



experiment has reported a positive signal ||26l . They used Nal, in which both nuclei 
have spin, one with an unpaired proton and the other with an unpaired neutron. 
The interpretation of their signal in terms of a WIMP with scalar interactions was 
ruled out by null results (at the time) from CDMS. An interpretation of their signal 
in terms of a spin-dependent WlMP-neutron interaction was ruled out by the null 
search in their Xe detector ll27l . While the interpretation in terms of spin-dependent 
WlMP-proton scattering was consistent with null results from other direct searches 
ll27l . it was ruled out by null searches for energetic neutrinos from the Sun (see 
Fig.[8]l, as we explain below. The interpretation in terms of spin-dependent scattering 
is now also ruled out directly by null results from the COUPP experiment 1281 . 



16 Guido D'Amico, Marc Kamionkowski and Kris Sigurdson 

4.3 Energetic v's from the Sun 



The escape velocity at the surface of the Sun is ^ 600km/s, while at the cen- 
ter it is Vc ~ 1300km/s. If in passing through the Sun, a WIMP from the Galac- 
tic halo scatters from a nucleus (most likely a proton) therein to a velocity less 
than the escape velocity, then it is gravitationally trapped within the Sun. As the 
gravitationally-trapped WIMP passes through the Sun subsequently, it loses en- 
ergy in additional nuclear scatters and thus settles to the center of the Sun. In this 
way, the number of WIMPs in the center of the Sun is enhanced. These WIMPs 
can then annihilate to standard model particles, through the same early-Universe 
processes that set their relic abundance 1291 . Decays of the annihilation prod- 
ucts (e.g., W^W^ ,Z^Z^ ,tt,bt>,cc, . . .) to neutrinos will produce energetic 
neutrinos that can escape from the center of the Sun. The neutrino energies are 
Ey ~ [(1 /3) — (1 /2)]M^ ~ lOOGeV and so cannot be confused with ordinary so- 
lar neutrinos, which have energies ~MeV. At night, these neutrinos will move up 
through the Earth. If the neutrino produces a muon through a charged-current inter- 
action in the rock below a neutrino telescope (e.g., super-Kamiokande, AMANDA, 
or IceCube), the muon may be seen. In this way, one can search for these WIMP- 
annihilation-induced neutrinos from the Sun. 



4.4 Cosmic rays from DM annihilation 

In the Galactic halo, one expects the annihilation processes XX ^ * ^iPP:"/"/'^ 

detection of these products can be a signal of the presence of dark matter. 

Exercise 9. Show that the annihilation process XX ^ e^e^ is suppressed forMajo- 
rana WIMPs as the relative velocity v ^ 0. 



4.4.1 Positrons 

Because of Galactic magnetic fields, cosmic -ray positrons and antiprotons do not 
propagate in straight lines and will thus appear to us as a diffuse background. Con- 
tinuum e+'s from WIMP annihilation are difficult to separate from ordinary cosmic- 
ray positrons. It has been argued that indirect processes, such as the annihilation 
XX ^ W^W^ e+ve^V | |30| , will produce a distinctive bump in the positron 
spectrum at energies Eg < (direct annihilation of Majorana WIMPs to electron- 
positron pairs is suppressed at Galactic relative velocities), as illustrated in Fig. |9] 
and there has been tremendous excitement recently with the reported detection by 
the PAMELA experiment of such a bump |f3ll. However, it may be that nearby 
pulsars can also produce a bump in the positron spectrum [1321 . and more recent re- 
sults from the Fermi Telescope 1331 call the PAMELA result into questions. It will 
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Fig. 8: The shaded region shows the parameter space (in WIMP mass versus SD WlMP-proton 
cross section) implied by the DAMA annual modulation for a WIMP with exclusively SD interac- 
tions with protons and no interaction with neutrons. The solid curve indicates the upper bound to 
the SD WlMP-proton cross section from null searches for neutrino-induced upward muons from 
the Sun; thus the cross hatched region is excluded 1 27 1 . 



thus be important to understand the possible pulsar signal, as well as the data, more 
carefully before the PAMELA excess can be attributed to WIMP annihilation. 



4.4.2 Antiprotons 

Likewise, it has also been argued that low-energy antiprotons from WIMP annihi- 
lation can be distinguished, through their energy spectrum, from the more prosaic 
cosmic -ray antiprotons produced by cosmic-ray spallation. Antiprotons can be pro- 
duced by the decay of the standard WIMP-annihilation products, and the energy 
spectrum of such antiprotons is relatively flat at low energies. On the other hand, the 
energy spectrum of low-energy cosmic -ray antiprotons due to cosmic -ray spallation 
decreases at energies E <Ge V. This is because the process p + pisM ^ P+P + P + P 
has an energy threshold, in the center of mass, of Eqm > 4-mp- This requires the 
primary cosmic -ray momentum to be very high. Production of an antiproton with 
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Ep <GeV therefore requires that the antiproton be ejected with momentum opposite 
to that of the initial cosmic-ray proton, and the phase-space for this ejection is small. 



4.4.3 Gamma rays 



A final channel to observe WIMP annihilation is via gamma rays from WIMP an- 
nihilation. Direct annihilation of WIMPs to two photons, XX ^ YY^ through loop 
diagrams such as those shown in Fig. [TOl produce monoenergetic photons, with 
energies equal to the WIMP mass. For v ^ lO^^c, the photon energies would be 
Ey = E-j^ (^1 ± 10^'^), and one would see a narrow y-ray line with Av/v ^ 10^^, 
superposed on a continuum spectrum produced by astrophysical processes; such a 
line would be difficult to mimic with traditional astrophysical sources. Decays of 
WIMP-annihilation products also produce a continuum spectrum of gamma rays at 
lower energies. 




POSITRON ENERGY (GeV) 



Fig. 9: The positron fraction, as a function of electron-positron energy, from annihilation of a 
120 GeV neutralino WIMP to gauge bosons. The different curves are for different cosmic-ray- 
propagation models, and in both cases, the annihilation rate has been boosted by a factor of ten 
relative to the canonical (smooth-halo) value. From Ref. 1301 . 
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Fig. 10: Example of a Feynman diagram for annihilation of two neutralinos to two photons through 
a quark-squark loop. 



The other advantage of gamma rays is that they propagate in straight Hnes. This 
opens up the possibiUty to distinguish gamma rays from WIMP annihilation from 
those from traditional sources through directionality — there should be a higher flux 
of WIMP-annihilation photons from places where WIMPs are abundant; e.g., the 
Galactic center Another possibility is dwarf galaxies, which represent regions of 
high dark-matter density in the Milky Way halo. In general, the y-ray flux (the num- 
ber of photons per unit time-area-solid-angle) is given by 

dn AkMI Jo ^ ' ^ ' 

where the integral is taken along a given line of sight, I is the distance along that line 
of sight, and p [l) is the dark-matter density at that distance. (Note that if p (r) °^ r^' 
with Galactocentric radius r, as in an NFW profile, the intensity formally diverges, 
but the flux form any finite angular window around r = is finite.) 

Exercise 10. Estimate the y-ray flux from WIMP annihilation, for a given annihila- 
tion cross section (times relative velocity) (<7v)ann, in an angular window of radius 
^ 5 degrees around the Galactic center. Estimate a characteristic (<7v) for WIMPs 
and evaluate your result for the gamma-ray flux for that value. How does it compare, 
in order of magnitude, with the sensitivity of the Fermi Gamma Ray Telescope? 



4.4.4 Galactic Substructure and Boost Factors 

The rate for annihilation, per unit volume, at any point in the Galactic halo is pro- 
portional to p^, the square of the density at that point. The total annihilation rate in 
the halo, or in some finite volume of the halo, is then proportional to / dV p^, the 
integral, over that volume, of the density squared. In the canonical model, the halo 
density is presumed to vary smoothly with position in the Galaxy with some density 
profile; e.g., the isothermal profile in Eq. (|2|i. 

However, a Galactic halo forms as part of a recent stage in a sequence of hi- 
erarchical structure formation. In this scenario, small objects undergo gravitational 
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collapse first; they then merge to form more massive objects, which then merge to 
form even more massive objects, etc. If some of these substructures remain partially 
intact as they merge into more massive halos, then any given halo (in particular, 
the Milky Way halo) may have a clumpy distribution of dark matter This is in fact 
seen in simulations. What this implies is that the annihilation rate in the halo may 
be enhanced by a "boost factor" B °^ (p^) /(p)^, where the averages are over vol- 
ume in the halo 1341 . It may be possible to see angular variations in the gamma-ray 
signal from WIMP annihilation, due to this substructure 13511361 . It has even been 
suggested that proper motions of nearby substructures may be visible fyj], although 
Ref. [381 disputed this claim. 

As we will see below, the first gravitationally-collapsed objects in WIMP models 
have masses in the range 10^^ — 100 Earth masses f39l. These objects may have 
densities several hundred times those of the mean halo density today. If so, and if 
these Earth-mass substructures survive intact through all subsequent generations of 
structure formation, then the boost factor B may be as large as several hundred, 
implying much larger cosmic -ray fluxes than the canonical model predicts. 

Such large boost factors are, however, unlikely. Simulations of recent genera- 
tions in the structure-formation hierarchy show that while the tightly bound inner 
parts of halos may survive during merging, the outer parts are stripped. Ref. l40l 
developed an analytic model, parametrized in terms of a halo-survival fraction, to 
describe the (nearly) scale-invariant process of hierarchical clustering. This model 
then provided the boost factor B in terms of that survival fraction. By comparing 
the results (cf.. Fig. fTTI) of the analytic model for the local halo-density probabil- 
ity distribution function with subsequent measurements of the same distribution in 
simulations (Fig. 1 in Ref. BTl ). one infers a small halo-survival fraction. The ana- 
lytic model of Ref. BOl then suggests for this survival fraction no more than a small 
boost factor, B <few. 



5 Variations and additions 

What we have described so far may be referred to as the minimal-WIMP scenario. 
In this scenario, the dark matter is a thermal relic with electroweak-scale cross sec- 
tions. It is neutral and scatters from baryons with cross sections ^ 10^'^'^cm^ (to 
within a few orders of magnitude). It has no astrophysical consequences in the post- 
freezeout Universe beyond its gravitational effects. However, the recent literature is 
filled with a large number of astrophysical anomalies for which explanations have 
been proposed in terms of non-minimal WIMPs, WIMPs endowed with extra in- 
teractions or properties. This is a vast literature, far too large to review here. We 
therefore provide here only a brief sampling, focusing primarily on those that we 
have worked on. 
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5.1 Enhanced relic abundance 

The calculation above of the freezeout abundance is the standard one in which it is 
assumed that the Universe is radiation-dominated at 7/ ^ 10— lOOGeV. However, 
we have no empirical constraints to the expansion rate before big bang nucleosyn- 
thesis, which happens later, at Tbbn ^ 1 MeV. 

One can imagine other scenarios in which the WIMP abundance changes. For 
instance, suppose the pre-BBN Universe is filled with some exotic matter which has 
a stiff equation of state, ps = Ps- This results in a scaling of the energy density of 
this stuff Ps °^ a^^ with scale factor a ll42ll . Such an equation of state may arise if the 
energy density is dominated by the kinetic energy of some scalar field. The equation 
of motion of a scalar field with a flat potential is 



(p + 3Hf = 



In^ cx: —3 In a 



(16) 
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which means that 

P = i(p2c<fl-6. (17) 

A stiff equation of state, or something that behaves effectively like it, may also arise, 
for example, in scalar-tensor theories of gravity or if there is anisotropic expansion 
in the early Universe. 

Big-bang nucleosynthesis constrains the energy density of some new compo- 
nent of matter at a temperature T ~MeV to be (pe/Py) ^ 0.1 (T/MeV)". Since 
Ps/ Pad °^ the expansion rate with this new stiff matter will at earlier times be 
H{T) < Hst{T) (r/MeV), where H^[{T) is the standard expansion rate. Neglect- 
ing the logarithmic dependence of the freezeout temperature Tf o= ln[Hp(,ny\ on the 
expansion, the WIMP abundance with this new exotic matter will be 

n_x^lL^LIL<('ll) (J-)^('^) (^). (18) 
«y «y (TV ny av ^ \ny J \MeV J \«y/^., \ MeV / 

Thus, for example, the relic abundance of an ~ 1 50 Ge V WIMP can be increased 
by as much as ~ 10'* in this way t42u43i . 

Exercise 11. Show that anisotropic expansion gives rise to a Friedmann equation 
that looks like that for a Universe with a new component of matter with p °^ a^^. 
To do so, consider a Universe with metric ds^ = dt^ — [aj{t)]^dx^ — [ay{t)]^dy-^ — 
[a-{t)]^dz^, with ax{t), ay{t), and az{t) different, and then derive the Friedmann 
equation for a Universe filled with homogeneous matter of density p . 



5.2 Kinetic decoupling 

There are two different kinds of equilibrium for WIMPs in the primordial bath. One 
is chemical equilibrium, which is maintained by the reactions 

the other is kinetic equilibrium, maintained by the scattering 

Xf^Xf- 

The first reaction freezes out before the second, since «/ S> n^, where / is 
any kind of light degree of freedom. However, <y{vx ^ vx) °= since the v's 
are Yukawa coupled, and (J{yX ^ YX) °^ Ey since the photons are coupled by 
e^ivpak^k^ £^ Ell- This means that r[xf <-> xf) drops rapidly and so kinetic 
freezeout happens not too much later than chemical freezeout. 

Detailed calculations of the kinetic-decoupling temperature 7].^ show that 
varies over 6 orders of magnitude in scans of the SUSY and UED parameter 
spaces 1391 . During the time particles are chemically but not kinetically decou- 



Dark Matter Astrophysics 



23 



pled, they have the same temperature of the thermal bath, which scales as oc a \ 
and after that, T-^ ~ p^^jlM-^ °^ a^^. So, density perturbations 5px/Px are sup- 
pressed on Ap/,„ ~ while the WIMPs are kinetically coupled. The cutoff in 
the power spectrum P{k) is at physical wavenumber = H{Ti^^j), so if 7].j de- 
creases, also kc decreases. We expect power suppressed at mass scales M < Mc, 
where Mc ~ lO""* - lO^M^ is the mass enclosed in the horizon at T^^, as shown in 

Fig.iniiiiii. 

Exercise 12. Derive the mass M^d enclosed within the horizon at a temperature T/^. 
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Fig. 12: The wavenumber and mass scale at which the primordial power spectrum is cut off due to 
kinetic decoupling of WIMPs in supersymmetric and UED models for WIMPs. From Ref. 1391 . 
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5.3 Particle Decay and Suppression of Small Scale Power 

It might be the case that dark matter is produced by the decay of a metastable particle 
that was once in kinetic equihbrium with the thermal bath. For instance, although 
the dark matter cannot be a charged particle it might be produced by the decay of 
a charged particle. The growth of perturbation modes that enter the horizon prior to 
the decay of the charged particle will be suppressed relative to the standard case due 
to the coupling to the thermal bath: growth of charged-particle density perturbations 
is suppressed since charged particles cannot move through the baryon-photon fluid. 
If one has x'^ ^ z'' + with T ^ 3. Syr (z ^ 10^), then the matter power spectrum 
P{k) is suppressed onk> Mpc^' ||45]| . while for shorter lifetimes structure will be 
suppressed for larger k (smaller length scales). Models exhibiting charged-particle 
decay can be found in the parameter space of standard or minimal extensions of 
canonical WIMP (e.g., supersymmetric) scenarios B6l . While limits on energy in- 
jection and the formation of exotic bound states in big bang nucleosynthesis (BBN) 
constrain the fraction of the Universe bound up in charged particles BtI the suppres- 
sion of power due to particle decay in the Universe remains a potentially observable 
effect of metastable particles. It is possible the metastable particle might remain in 
kinetic equilibrium via another interaction, or even if the particle is out of kinetic 
equilibrium the energy released in the decay process may impart the dark-matter 
particle with a velocity high enough to erase small-scale structure via free stream- 
ing iU. 

Future measurements of high-redshift cosmic 21 -cm fluctuations may provide 
a direct probe of modifications to the small-scale dark-matter power spectrum and 
other aspects of fundamental physics (see, e.g., 14611491 ). 

Exercise 13. Derive the comoving wavenumber k that enters the horizon at the time 
a particle of lifetime T decays. 



5.4 Dipole dark matter 

While dark matter cannot be a charged particle it may (via higher order interactions) 
be endowed with an electric or magnetic dipole moment interactions of the form 

Ldipoie - Xi(y^iv i^ij + Ys^ij) XjF'"", (19) 

Here, diagonal interaction terms (; = j) are the magnetic (/i) or electric (^) dipole 
moments of a particle while off-diagonal terms (/ ^ j) are referred to as tran- 
sition moments between the lightest WIMP state / and another, slightly heavier, 
WIMP state j. Such a dipole coupling to photons alters the evolution of dark-matter 
density perturbations and CMB anisotropies ll22ll . although the strongest constraints 
to dipole moments comes from precision tests of the Standard Model for WIMP 
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Fig. 13: Constraints to the dipole-mass parameter space for dark matter with an electric or magnetic 
dipole. From Ref. 1221 . 



masses My, < 10 GeV and direct-detection experiments for My, > 10 GeV 
see Fig.[T3]for the full constraints. 

It may be possible to explain the results of the DAMA experiment using low- 
mass dipolar dark matter with a transition moment ll50l . It may also be possible to 
look for the effects of a transition dipole moment in the absorption of high energy 
photons from distant sources IT91. 



Exercise 14. Calculate the cross section for elastic scattering of a particle with an 
electric dipole moment of magnitude d from a nucleus with charge Ze. 



5.5 Gravitational constraints 



It is generally assumed that while dark matter may involve new physics, the gravi- 
tational interactions of the dark matter are standard. In other words, it is generally 
assumed that the gravitational force between two DM particles and between a dark- 
matter particle and a baryon is the same as that between two baryons. More pre- 
cisely, the Newtonian gravitational force law between baryons that has been tested 
in the laboratory and in the Solar System reads Ft,^i,2 = Gm\m2/d~. We then usually 
assume that the force between baryons and DM is Fhd = Gmi,mj / d^, and also that 
the gravitational DM-DM force law is /Vic/i = Gdindiinij^/d^ with 0^ = G. How- 
ever, there is no empirical evidence that this is true at more than the order-unity level 
ifSTl . and it has even been postulated that = 2G in order to account for the void 
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abundance ||52l . A similar behavior (an increase in the DM-DM force law) could 
also arise if there were a new long-range interaction mediated by a nearly massless 
scalar field (p with Yukawa interactions (pxjfitf with the DM field The difficulty 
in providing empirical constraints to this model is that measurements (e.g., gravita- 
tional lensing or stellar/galactic dynamics) of the dark-matter distribution determine 
only the gravitational potential <P due to the dark-matter distribution, represented by 
some density Pdif)^ obtained through the Poisson equation V^^> = A-KGpd- How- 
ever, the same <P can be obtained by replacing Pd ^ {l/2)pd if we simultaneously 
replace G — > 2G. 

It turns out, though, that this exotic interaction can be constrained by looking 
at substructures in the Milky Way halo ll53l 1541 . The Sagittarius dwarf galaxy, is 
dark-matter dominated, and it follows an elongated orbit around the Milky Way. 
When the dwarf reaches its point of closest approach to the Milky Way, the tidal 
forces it experiences in the Milky Way potential are largest. Stars are then stripped 
from the innermost and furthermost edge of the dwarf. Those from the innermost 
parts move at slightly larger velocities in the Galactic halo and at slightly smaller 
Galactocentric radii; they thus subsequently run ahead of the Sagittarius dwarf and 
form the leading tidal tail of the Sagittarius dwarf that is observed. Conversely, those 
stripped from the outer edge subsequently lag behind forming the trailing tidal tail 
that is observed. Observationally, the leading and trailing tails have roughly the same 
brightness, as expected. Suppose now that the DM-DM force law were modified to 
Gd = fG with / > 1. The dark-matter halo of the Sagittarius dwarf would then be 
accelerated toward the Milky Way center more strongly than the stellar component 
of the Sagittarius dwarf. The stellar component would then slosh to the furthermost 
edge. Then, when the dwarf reaches its point of closest approach to the Milky Way, 
stars are still stripped from the outer edge, forming a trailing tail. However, there 
are now no stars in the innermost edge to form the leading tail. The evacuation of 
stars from the leading tail is inconsistent with observations, and this leads, with 
detailed calculations, to a bound G^/ = G(l ±0.1) to Newton's constant for DM- 
DM interactions. In other words, dark matter and ordinary matter fall the same way, 
to within 10%, in a gravitational potential well. 

While Ref. Il55l has more recently claimed to run a simulation of the tidal tails of 
the Sagittarius dwarf consistent with Gd = 2G, Ref. B56l has argued that the initial 
conditions for that simulation are self-inconsistent. Refs. Il57l l58l argue that a new 
long-range DM-DM force law implies, under fairly general conditions, a weaker 
long-range DM-baryon force law, and they discuss and compare possible tests of 
such a scenario. 



5.6 Electromagnetic-like interactions for dark matter? 

Another possibility is that dark matter experiences long-range electromagnetic-like 
forces mediated by a dark massless photon that couples only to gravity. Of course, if 
the fine-structure constant associated with this dark U{1) symmetry is too large. 
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then long-range dark forces will induce the dark matter to be effectively collisional. 
This constrains < 0.005 {M^/TeW)^/^ ll59l . Far more restrictive constraints may 
arise from the development of plasma instabilities that may arise if there are (dark) 
positively and negatively charged dark-matter particles, but precise calculations of 
these effects remain to be done. See Refs. |f59l l60l for more discussion of these 
models. 

Exercise 15. Estimate the relic abundance of a dark-matter particle with dark charge 
ad assuming that it annihilates to dark-photon pats and assuming that the dark 
sector has the same temperature as the rest of the primordial plasma. 



6 Some other particle dark-matter candidates 

WIMP models are interesting for a number of reasons: (1) The correct relic density 
arises naturally if there is new physics at the electroweak scale; (2) there are good 
prospects for detection of these particles, if they are indeed the dark matter; and (3) 
there is synergy with the goals of accelerator searches (especially at the LHC) for 
new electroweak-scale physics. 

Still, there are a large number of other particle candidates for dark matter. Here 
we discuss two, the sterile neutrino and the axion, which may also arise in extensions 
of the standard model and for which there are clear paths toward detection if they 
make up the dark matter 



6.1 Sterile Neutrinos 

A convenient mechanism to introduce neutrino masses and explain their smallness 
by a minimal extension of the Standard Model is to add 3 right-handed neutrinos 
which are singlets under the SM gauge group. The mass matrix is taken to be of the 
form (for simplicity we consider only one family). 



where the Vl and Vr are left-handed and right-handed (weakly-interacting and ster- 
ile, respectively) fields. 

In the see-saw mechanism, the Dirac mass is assumed to be tiny compared with 
the Majorana mass: i.e.. Mo ^ M. The mass eigenstates then have masses Mi ~ 
M^/M <^ M, and M2 — M. For our purposes, it is advantageous to map the two- 
dimensional Mjj-M parameter space onto the M^-O parameter space, where Ms is 
the mass of the sterile (heavier) neutrino and 9 is the mixing angle between the two 



Vl Vr 




(20) 
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States. The active and sterile mass eigenstates can then be written 

|v„) =cos0|Vi)+sin0|Vfl), (21) 
|v,) = -sin0|vL)+cos0|Vff), (22) 

where 9 ^Mo/M. 

Sterile neutrinos can be produced in the early Universe and have both (1) a life- 
time longer than the age of the Universe and (2) a cosmological density £2^ ^ 0.2 if 
the sterile-neutrino mass is in the ~keV regime IMI . 

The main decay mode of the sterile neutrino is then Vg VVV, through the 
exchange of a boson, as shown in Fig. [14] The decay rate and lifetime are 




Fig. 14: Main decay channel for sterile neutrinos. 

If the sterile neutrinos constitute the dark matter, then it must be that Ts ^ 
lO^^sec, which is possible if Ms ^ (9(l)keV. This mass cannot however be too 
small, because of the Gunn-Tremaine limit from dwarf-spheroidal galaxies, which 
is Ms > 0.3 keV. A stronger constraint to the model comes from the X-ray emis- 
sion in the radiative decay Vs vy, through the diagram in Fig. [15] This produces 
an x-ray line that can be sought in the spectrum of, e.g., a galaxy cluster. While 
null searches for such lines (and from the diffuse cosmic x-ray background) pro- 
vide ll62l |63]| stringent constraints to the model, there are still some regions in the 
Ms-9 parameter space that remain consistent with current constraints. This region 
may be probed, however, with future more sensitive x-ray searches. One interesting 
extended application of sterile neutrino dark matter was its use as a potential mech- 
anism for generating momentum-anisotropy during supernova to drive pulsar kicks 
li64ll . See, for instance, Refs. li65]|66l . for the current status of sterile neutrino dark 
matter. 
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Fig. 15: Loop diagram for the decay v, — > vy. 

6.2 Axions 

Axions arise in the Peccei-Quinn (PQ) solution to the strong-CP problem 1671 . A 
global U{l)pQ symmetry is spontaneously broken at a scale /„, and the CP- violating 
phase 9 in the QCD Lagrangian becomes a dynamical field with a flat potential. 
At temperatures below the QCD phase transition, nonperturbative quantum effects 
break explicitly the symmetry and produce a non-flat potential that is minimized at 
0^0. The axion is the pseudo-Nambu-Goldstone boson of this near-global sym- 
metry, the particle associated with excitations about the minimum at = 0. The 
axion mass is rria — eV (lO^GeV/Zo), and its coupling to ordinary matter is °^ f^^ ■ 

The Peccei-Quinn solution works equally well for any value of /„. However, a 
variety of astrophysical observations and laboratory experiments constrain the axion 
mass to be nia ^ 10^^ eV. Smaller masses would lead to an unacceptably large cos- 
mological abundance. Larger masses are ruled out by a combination of constraints 
from supernova 1987 A, globular clusters, laboratory experiments, and a search for 
two-photon decays of relic axions. 

Curiously enough, if the axion mass is in the relatively small viable range, the 
relic density is Q.^, ^ \, and so the axion may account for the halo dark matter 
Such axions would be produced with zero momentum by a misalignment mecha- 
nism in the early Universe and therefore act as cold dark matter. During the process 
of galaxy formation, these axions would fall into the Galactic potential well and 
would therefore be present in our halo with a velocity dispersion near 270 km sec^ ' . 

It has been noted that quantum gravity is generically expected to violate global 
symmetries, and unless these Planck-scale effects can be suppressed by a huge fac- 
tor, the Peccei-Quinn mechanism may be invalidated 1681 . Of course, we have at 
this point no predictive theory of quantum gravity, and several mechanisms for for- 
bidding these global-symmetry violating terms have been proposed |l69| . Therefore, 
discovery of an axion might provide much needed clues to the nature of Planck-scale 
physics. 

There is a very weak coupling of an axion to photons through the triangle 
anomaly, a coupling mediated by the exchange of virtual quarks and leptons. The ax- 
ion can therefore decay to two photons, but the lifetime is Ta^yy ^ 105°s(m„/10-5eV) 
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which is huge compared to the lifetime of the Universe and therefore unobservable. 
However, the ayy term in the Lagrangian is ^uyy °^ aE ■ B where E and B are the 
electric and magnetic field strengths. Therefore, if one immerses a resonant cavity in 
a strong magnetic field. Galactic axions that pass through the detector may be con- 
verted to fundamental excitations of the cavity, and these may be observable ifTOll . 
Such an experiment is currently underway ItTI and has already begun to probe part 
of the cosmologically interesting parameter space (see the Figure in Ref. 172]), and 
it should cover most of the interesting region parameter space in the next few years. 

Axions, or other light pseudoscalar particles, may show up astrophy sic ally or ex- 
perimentally in other ways. For example, the PVLAS Collaboration 1731 reported 
the observation of an anomalously large rotation of the linear polarization of a laser 
when passed through a strong magnetic field. Such a rotation is expected in quantum 
electrodynamics, but the magnitude they reported was in excess of this expectation. 
One possible explanation is a coupling of the pseudoscalar FF of electromagnetism 
to a low-mass axion-like pseudoscalar field. The region of the mass-coupling param- 
eter space implied by this experiment violates limits for axions from astrophysical 
constraints, but there may be nonminimal models that can accommodate those con- 
straints. Ref. ||74| reviews the theoretical interpretation and shows how the interac- 
tions of axions and other axion-like particles may be tested with x-ray re-appearance 
experiments. While the original PVLAS results have now been called into question 
Ref. Il75l . variations of the model may still be worth investigating. 



7 Conclusions 

Here we have reviewed briefly the basic astrophysical evidence for dark matter, 
some simple astrophysical constraints to its physical properties, and the canoni- 
cal WIMP model for dark matter. We then discussed a number of variations of the 
canonical model, as well as some alternative particle dark-matter candidates. Still, 
we have only scratched the surface here, surveying only a small fraction of the pos- 
sibilities for non-minimal dark matter. Readers who are interested in learning more 
are encouraged to browse the recent literature, where they will find an almost end- 
less flow of interesting possibilities for dark matter, beyond those we have reviewed 
here. 
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